Plot the consensus tree
tree_global = read.tree('RAxML_bipartitions.result_FIN4_raw_rooted_wBoots_4098mam1out_OK.newick')
n_tips = length(tree_global$tip.label)
root_node = tree_global$edge[(!tree_global$edge[,1] %in% tree_global$edge[,2]),1] %>% unique()
root_age = max( dist.nodes(tree_global)[,root_node] )
label_offset = 0.05 * root_age
tree_width = 1.5 * root_age
ggtree_global = ggtree(tree_global, lwd = 0.1) +
geom_tiplab( fontface = "italic", offset = label_offset, size = 2) +
xlim(0, tree_width)
ggtree_global
